Amplification of endpoint structure for new particle mass measurement 

at the LHC 
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We introduce a new collider variable, Mcti, named as constransverse mass. It is a mixture 
of 'stransverse mass(MT2)' and 'contransverse mass(AfcT)' variables, where the usual endpoint 
structure of Mt2 distribution can be amplified in the Mct2 basis by large Jacobian factor which 
is controlled by trial missing particle mass. Thus the Mcti projection of events increases our 
observability to measure several important endpoints from new particle decays, which are usually 
expected to be buried by irreducible backgrounds with various systematic uncertainties at the LHC. 
In this paper we explain the phenomenology of endpoint amplification in Mcti projection, and 
describe how one may employ this variable to measure several meaningful mass constraints of new 
particles. 
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I. INTRODUCTION AND MOTIVATION 

The large hadron collider(LHC) will explore the TeV 
scale soon searching for new physics beyond the Standard 
Model(SM) 0,0, and various new physics models are 
waiting to be tested. 

A well-known expectation about the new physics phe- 
nomenology at the LHC is that a Zi parity conservation 
could be a common feature in some new physics mod- 
els where the lightest new particle(LNP) are very stable, 
providing weakly interacting massive particle(WIMP) 
dark matter candidates. Supersymmetry(SUSY) with Im- 
parity [H, Little Higgs models with T- parity or Uni- 
versal Extra Dimension models with the Kaluza-Klein 
parity [f| are the examples of such models. In those 
models, a pair of weakly interacting stable LNPs are 
missing in the detector leaving rich missing transverse 
energy signals. In general, the existence of multi- 
ple missing particles makes an event reconstruction very 
hard at the LHC where we will also suffer from the par- 
tonic center of momentum(CM) frame unambiguity and 
complex new event topologies. Under this circumstance, 
the mass measurement of new particles is not an easy 
task at the LHC, and many important previous studies 
on the subject might be useful [6H24|. 

In particular, when a new event cannot be recon- 
structed, measuring the endpoints of event projections 
onto various observables can be a good way to obtain 
some mass constraints of the new particles involved. It 
is because the endpoints correspond to the kinematic 
boundaries of the allowed phase space of the event, which 
arc mainly described by the related particle masses. 
Some examples of such kind of the methods include kine- 
matic endpoint methods using invariant masses of visible 
particles [7j , and mass measurements using the endpoints 
of various transverse mass variables [8H15| . The existence 
of such endpoints or some cusp points in a distribution 
can be re-analyzed by surveying singularity structures of 



the allowed phase space of the events [16| . 

For these methods using endpoints, the precise mea- 
surement of the endpoint must be the most important 
factor for it to be reliable. However, identifying the cor- 
rect endpoint is not an easy task in the real situation 
since in general there exist complex systematic uncertain- 
ties. Once we precisely understand the dynamics of the 
signal/background processes and its detector responses, 
then the mass parameters can be obtained via the least 
likelihood method using templates as in Ref. @. This 
method is reliable because the Standard Modcl(SM) has 
been well-understood, explaining most of the phenom- 
ena in the high energy experiments up to now. However, 
when it comes to the LHC, things will be changed. If 
there exist new physics beyond the SM at the TeV scale, 
it is quite obvious that we will suffer from huge systematic 
uncertainties in measuring the masses or new model pa- 
rameters. For example, about the signal with jets, heavy 
jet combinatoric backgrounds should be understood in 
the situation where the new physics events can have very 
complex event topologies with hard QCD effects on the 
beyond the SM(BSM) signals. 

In this regard, extracting the meaningful endpoints via 
some simplified functional fitting have been performed 
near a proper and plausible endpoint region of the dis- 
tribution [7h15|. In that region, a Gaussian smeared 
segmented linear functional fitting for signal and back- 
grounds is likely to give a sufficiently good description 
in many cases. However, it is true that one would still 
suffer from large systematic uncertainties in these seg- 
mented straight-line fits to extract the endpoint position. 
In particular, it is more severe when the signal endpoint 
structure is faint with feet or tails (ignoring any smearing 
effects from finite total decay width or detector resolu- 
tions), leaving small number of events near the endpoints, 
not constructing a sharp drop [l7| . In many cases, the 
endpoint fitting has quite a large uncertainty with irre- 
ducible backgrounds. 
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In this paper, we assume a situation where one suffers 
from heavy systematic uncertainties in identifying the 
cndpoints, and define the meaningful signal endpoint as a 
breakpoint (BP) in the distribution up to smearing effects. 
We describe a way to obtain some mass constraints more 
precisely, which have been obtained from the endpoints 
of Mt2 distributions. To do so, we introduce a new vari- 
able, Mcti, named as 'constrans verse mass'. The Mcti 
projection of the events is found to have interesting Ja- 
cobian factor with respect to Mj- 2 distribution, which 
amplifies the endpoint structures of My 2 . 

In Sec. II, we start with the definition of Mct which is 
the basic ingredient for the Mcti study, and describe the 
phenomenology of the endpoint structure amplification, 
comparing it with the Mr distribution. In Sec. IIII| the 
Mcti variable is introduced. In Sec lIVI it is described 
how one may employ the Mcti projection for a precise 
measurement of the signal endpoint, usually buried in 
the irreducible jet combinatoric backgrounds. Sec|V]is a 
conclusion of our study. 

II. PROPERTIES OF M C t 

Let us consider a mother particle Y decaying to the 
visible and missing particles (Y — > v(p) + X(k)), where p 
and k are their corresponding 4 momenta. Then, one can 
reconstruct transverse mass(MT) Q which is bounded 
from above by the invariant mass(M) of Y, 

M 2 = ml + m\ + 2{E V E X -p-k) 

= ml + m\ + 2(e v e x cosher?) - p T ■ k T ) (1) 
> Mf 

where e v = \J ml + \pt | and ex = \J m 2 x + | A;t | 2 are 
transverse energies of the visible and missing particles 
with their transverse momenta, pr and kr, A 77 is the 
rapidity difference between v and X, and Mf is defined 
by 

Mf = ml +m 2 x +2(e v e x - pr • fer). (2) 

Furthermore, we can also define Mc and Mct in terms 
of v and X pl], 

M c = m 2 v + m 2 x + 2A, 

M CT = m 2 v + m 2 x + 2 A T , 

~ (3) 
A = E V E X +p-k, 

A T = e v e x +Pt ■ k T , 

satisfying a similar inequality, 

M 2 > M 2 T . (4) 

An interesting property of Mc variable is that it is in- 
variant under the back-to-back boost of the two particles, 

p"->ASGdy, k^K{-p)k v , (5) 



where A^(/3) denotes the Lorentz transformation matrix 
for the boost parameter (3. The back-to- back boost in- 
variance(BBI) results from the BBI of "A-term", of 
which is the Euclidean dot products of p and k. The BBI 
of the Euclidean momenta product has been noticed in 
pll [20j with the form of At in Mti solutions, and uti- 
lized in Mct [HI using the visible transverse momenta. 
In the rest frame of Y, Mc and Mct read 

M 2 = ml + m 2 x + 2( e y x cosh(SV) - |p 0T | 2 ) 
M 2 T = m 2 v +m 2 x + 2(e° v e x -\poT\ 2 ) 

where e° x are the transverse energies of the daugh- 
ter particles, Hrf is the rapidity sum of v and X, and 
|po| = A 1 / 2 (my, m x , ml)/2mY, 1 which is the absolute 
momentum of daughter particles in the Y rest frame. 
The reason that we specify Mc and Mct in the rest 
frame of Y as in ©, is because Mc is not a frame inde- 
pendent quantity. From now on, we will use the value of 
Mc defined only in the rest frame of mother particle as 
in (0J), however we will not require such a restriction on 
M CT - 

The distribution of Mct is interesting. In particular, 
the range of the distribution can be drastically changed 
with respect to the input trial missing particle mass x- 
Let us assume that the true value of missing particle mass 
is not known. We can also construct Mt(x) an d Mct(x) 
using a trial invisible mass, %, by replacing mx to \ ln 
(|1|2|6|) . Then, as long as the Y is transversely at rest 
in the lab. frame, Mct(x) is bounded from above by 
Mc{x) as follows 

M 2 = X 2 + 2(\p \E x ~\p \ 2 ), 
M CT = x 2 + 2(| Po t|4-|pot| 2 ), (7) 

M% > M CT 

in which Mt(x) a l so satisfies a similar relation, 

M 2 = X 2 + 2{\p \E\ + | Po | 2 ) 
M 2 =x 2 + 2(| P oT|e° Y + | P oT| 2 ), (8) 
M 2 > Mf 

where Ex (ex) = \/x 2 + |Po (Pot)| 2 , and the visible 
masses are ignored in both cases. It is very well known 
that Mt{x = m x) has a nice endpoint, my, which is in- 
variant under the transverse boost of Y , even though the 
bulk distribution is shifted by Pt of Y. However, when X 
is not the true value (the maximum of Mt(x)), M( X ) is 
not invariant any more under the transverse motion of the 
Y 2 It is also true for the maximum of Mct(x m x)- 



1 X(x, y, z) = x 2 +y 2 + z 2 - 2(xy + yz + xz) 

2 The features of transversely boosted My (y) and Mj^Cx) end- 
points were surveyed by A. Barr et. al. in [1 j| and M. Burns et. 
al. in [TJ 
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In the following section, we will comment also on the 
feature of the shift of Mct{x) endpoint. 

Then, let us compare the distributions of Mqt(x) and 
My(x) when the mother particle is transversely at rest in 
the lab. frame. From Eqs. (|7I8[) . we note that Mct(x) 
and Mt(x) share the same minimum point at x but get 
different maxima. As x tends to zero, Mc{x) tends zero, 
and Mct{x) distribution collapses almost to a zero point, 
while Mt ranges from zero to 2|po|. On the other hand, 
when x is very large compared to |po|, their ranges and 
shapes become very similar because the effect of the sign 
flip in the squared momentum is negligible. The most 
impressive property of the Mct(x) projection is that 
if we use a proper value of x(= x/|Po|), which is not 
quite large, then the endpoint structure of Mct{x) dis- 
tribution can be enhanced by some large Jacobian factor 
Jmax(x) compared to the Mt{x) endpoint distribution. 
The Jacobian factor is given by 



da 



dM C r{x) 



(1(7 



dM T (x)' 



or 



J 



Mct(x) (ex + |Pot|) 2 



M t (x) ( e x 
whose behavior is 

M c (x) (Ex + \Po\ f _ 



I Pot | 



(9) 



M(x) (£x-|po|) 2 



1, 



J rt 



the endpoint region, 
the minimum region. 



Thus, the endpoint enhancement factor J max approaches 
oo when x(= x/|Po|) goes to zero, or becomes one if 
X is very large. As a result of the very different com- 
pression rate between the endpoint and minimum re- 
gion of Mt, most of the large Mt events are accumu- 
lated in the narrow maximal region of Mct distribution 
if x is not so large. Fig. Q] shows the distribution of 
Mct(x) and Mt(x) when Y is transversely at rest in 
the lab. frame with my = 100 GeV, m x = 10 GeV, 
and m v = GeV. The total decay width of Y, Ty, 
was assumed to be 1 GeV and spin angular correlation 
effect was ignored. In Fig. [Ua), we used a trial miss- 
ing particle mass x corresponding to 2|po|(= 99 GeV) 
where J max is 5.39. The black dashed line filled with 
green hatch corresponds to Mt(x) ~ X distribution and 
the black line is Mct(x) ~ X distribution. As explained, 
Mct and Mt share the same minimum at Xi but the 
endpoint of Mct shrinks into a lower region so that the 
expected endpoint is M c 1 ^ x (99) - 99 = 26.9 GeV while 
M^ ax (99) - 99 = 61.2 GeV, ignoring the width effect. 
The height of the peak in Mct distribution is amplified 
by the J m ax factor compared to the peak of Mt- Fig. 
[IJb) shows how Mct and Mt are varied with respect to 
the trial missing particle mass x- As shown in the figure, 
when the x is small, J ma x diverges(red: J ma x = 105), 
and Mct distribution collapses into a narrow peak re- 
gion. On the other hand, if x is large, then Mct and 
Mt become the same (black: J ~ 1). 
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FIG. 1: Enhancement at the endpoint. 



Now we have a collider variable Mct(x) which can 
drastically change and magnify its endpoint structure 
compared to that of Mt(x)- Basically, they provide the 
same physical constraint. Nevertheless, there is a big 
difference on the systematic uncertainties at these two 
end point regions. So, we move on to discuss how the 
amplification of the endpoint structure accompanies the 
propagation of statistical and systematic uncertainties. 

In principle, one can try to use the sharp amplified 
Mct{x) endpoint instead of the Mt(x) endpoint in order 
to measure the mass constraint |po| more precisely. How- 
ever, one must be sure about that a naive statistical un- 
certainty cannot be reduced since as SM™f x (x) decreases 
to -jl — 8M T ^ ax (x), the error propagation factor for |po| 
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increase by J m ax- On the other hands, things change 
when it comes to the systematic uncertainties. Actually, 
the endpoint structures we want to resolve are not that 
simple. If we attempt an accurate mass(endpoint) mea- 
surement with exact distribution shapes or templates, 
one have to understand the important dynamics of the 
new signal and backgrounds (surely also with the good 
understanding of the detector responses), which must 
include large systematic uncertainties. In this regard, 
model independent mass (endpoint) measurement will be 
important. 

An effective way to measure the endpoint is simply 
observing and pinpointing the breakpoint(BP) with sim- 
plified local fit function up to some smearing width in 
the distribution. In this respect, it is important to know 
the effective range near the BP for the simplified local fit 
functions to be reliable, because we suffer from a large 
uncertainty in the bulk distribution. However, it is not 
apparent when the signal endpoint structure is dim and 
faint with feet or long tail above some irreducible back- 
grounds so that the BP is quite ambiguous with a small 
slope discontinuity. [Wc will assume that the background 
distribution is not singular at the BP.] When we attempt 
to identify a BP as the position of a meaningful endpoint , 
a slope discontinuity between the lower and the upper re- 
gion of the BP is crucial factor for the BP measurement 
with less errors. The propagated uncertainty of the BP 
can be given as follows: 



Then 6%, 



1 ssys 



and if we take into account 



Aa 2 



(10) 



where a is an error which can be caused by various 
sources, and Aa is the slope difference between the two 
regions segmented by the BP. a may come from event 
statistics or some systematic uncertainties. The point 
we note here is that the larger the slope difference, the 
clearer the endpoint structure, and as a result it enables 
us to elaborate the fit function and to choose more effec- 
tive range of the fit function. In this regard, the elabo- 
ration of the fit scheme will reduce a sys , and as a result 
Sgp will be reduced. 

For the same reason an ambiguous faint endpoint in 
the Mt distribution can be measured more precisely in 
the Mct projection, and as a result the systematic er- 
ror in obtaining the mass constraint |po|, $p y S y can be 
reduced by 0(1/ J max)- The specific estimation of error 
suppression factor is based on the following two facts, 

(1) The Mct(x) projection amplifies the slope differ- 



ence Aa near the Mt(x) endpoint by J max , 



(11) 



(2) By the enhancement of the slope discontinuity, a sys 
in the Mct projection can be less than, or at least 
comparable to the Mt case with refined fit scheme 
near the enhanced endpoint 



the error propagation factor J ma x for |po|, we get S^ o s 



_J_Xsys 

4„ Po ■ 

Up to now we described that Mct(x) is well-defined 
variable using pt and trial missing particle mass \ m 
the rest frame of Y . The large amplification factor J max 
can be obtained simply by changing Xi an d there is no 
additional parameter needed to accentuate the dim BP of 
the Mt(x) distribution. One might have a good chance 
to reduce some systematic uncertainty in extracting the 
mass constraints with the Mct projection. 



III. CONSTRANSVERSE MASS, M C ti 

In this section wc introduce constransverse 
mass(Mc7T2) for new physics events with two miss- 
ing LNPs, which inherits the properties of Mct- Some 
differences of using Mct and Mct2 will be discussed 
also. 

When can we observe such an Mct endpoint in a 
real experiment? One thing we should remind is that 
the maximum of Mct(x) an d -^r(x), Mqj< x (x) and 
M™ ax (x), are not invariant under the mother particle's 
transverse motion when \ is n °t true value. One more 
shortcoming for M c ^ x (x) is that it is not also invariant 
under the transverse boost of the mother particle, even 
when x is true value. The reason for this is simply that 
Mc in ([3]) is not an invariant quantity under the boost 
of mother particle, while M(m x ) (the endpoint value of 
Mt(to x )) is an invariant one, as M™ ax has been uti- 
lized to measure the W-boson mass provided with already 
known neutrino mass 0, . However if we assume that 
we do not know the mass of the missing LNP, m x , while 
we just want to measure some mass constraints related to 
the decay process, then M^ x {x) has no disadvantages 
compared to M T n,ax (x)- As we will see in the examples 
of the following section, it can be much better to use 
M^ X ( X ) or Mg^f (x) for this purpose. Then, it is just 
required for us to choose our decay system of interest 
being at rest in the transverse direction. 

When the mother particle is boosted in the trans- 
verse direction x with momentum, 8t, the shift of the 
Mct{jti x ) endpoint is described by 



AM™" 



( m x) 



M c (m^ 



f(rn x ) a cos cj>, (a « 1) 



where f(m x ) 



4(Ex-l) 



Ex 



fsys 



(12) 



m|+2(B x -l)' ■"■> MPoF |Po 

a = St/tuy, and <j) is the azimuthal angle of visible par- 
ticle in the mother particle's rest frame. The shift is of 
order a in general, and one should take care of the use of 
the Mct(x) endpoint for the mother particle with sizable 
transverse momentum in the lab. frame. 5™ ax cut must 
be accompanied with a proper event selection. Surely, 
the St suppression cut might reduce the statistics of the 
signal event candidates. However, the cut can also play 
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a role for purifying the signatures from backgrounds in 
many cases. We will describe such an effect of the cut 
with a new physics example, using Mcti endpoint in the 
next section. 

Then, let us consider the production of a pair of identi- 
cal new physics particles at the LHC, in which each of two 
mother particles Ij decays to visible Vi and invisible par- 
ticles Xi, (p + p^ G(-S T ) + Yi{-^ Vi(pi) + Xi(ki)),i = 
1,2). Here G with its transverse momentum^— St) de- 
notes the other particles not from the decays, so that 
it provides the transverse momenta of Y\ + Y2 system as 
St- In this case, we find that it is also possible to ob- 
serve the same AfcT(x)-like endpoint by the construction 
of constransverse mass, Mct2(x) by 

M CT2 {x) = mm [max {AfH.M^.}] (13) 

M { cT la \x) 2 = X 2 + ™l, + 2{e Vi e Xi + p lT ■ k iT ), 

where x is a trial missing particle mass, and lfi T = 
— J2 PiT + St- Here all the particle momenta are denned 

in the lab. frame. M^ T (x) is the contransverse mass 
of each mother particle system as defined in (|3|) using a 
trial transverse momenta which satisfies the missing 
transverse momentum condition. This is the contra ver- 
sion of Cambridge strans verse mass (My 2) variable [nj. 
Please note that M^ T in ([T^| is defined for the mother 
particle pair where each mother particle can have sizable 
transverse momentum in the lab. frame, and it can make 
Mg) r ,nax (x) shifted as explained previously. However, an 
interesting point is that if St = 0, Y1+Y2 system is trans- 
versely at rest, and the M™ T lx {x) a lso has well-defined 
endpoint as Mc(x) m ©• This is analogous to the case 
that Mq^ x {x) becomes Mc(x) when the mother particle 
Y has no transverse momentum. To be analytically ex- 
act for that result, let us see the solution of Mct2(x)- If 
St = 0, the analytic Mct2(x) solution is given as follows, 

M 2 C T2(X) = x 2 -A t + ^AI + 2Atx\ (14) 
A T = |pit||P2t| + Pit • P2T, 

ignoring the mass of the visible particles. 3 The so- 
lution can be obtained just by using the visible mo- 
menta with the sign flip in solving the balanced equa- 
tion as in Refs. [ll|, [2fJ. From one can see that 
MffiZix) = M^ x (x) as the natural extension of M CT , 

M^UX) = M c {x) 

= X 2 +2(|p |^-|po| 2 ), 

in the limit of St = 0. This can be easily proved from 
Eq. (|14[) , assuming the situation where the two mother 



In the rest of this paper, we concentrate on the events with one 
masslcss visible particle in each decay chain. 



particles are at rest in the lab. frame and all the relic 
particles are in the transverse plane. Such a pair of rest 
mother particles can be guaranteed by the back-to-back 
boost invariance of Mcth(x) solution (| 14[) . It is interest- 
ing to note that M^^ix) 1S described by the po value 
in general, in spite of the nonzero transverse momen- 
tum for each of the mother particles. This is the same 
phenomenon that Mt2 inherits most of the important 
properties of M T , so that M^ x (x) = M^ ax (x) also. 
Physically Eq. (TT5|) can be satisfied by the events with 
kinematically identical decay chains, each of which has 
the maximal Mt or Mqt configuration. So in this case, 
Mt2 (or Mct2) can be effectively the same as the single 
Mt (01 Mct)- This is also true for the event contributing 
to the minimum. 

The other important inheritance of Mct2{x) from 
Mct(x) is that the Jacobian factor in © should be 
valid also between the Mct2(x) anc » Mt2(x)> a t least for 
the maximum and minimum regions of both variables as 
explained. In the maximum region, the Mct2(x) an d 
Mt2(x) solutions are perfectly described by Mc(x) m 
([T5"j) and M{x) in ©, respectively, while in the min- 
imum region, both Mct2(x) an d Mt2(x) become the 
same value, x gi vm g the Jacobian factor of 1. 

Fig. [5] shows how the previous analysis of Mt and 
Mct is extended to the Mt2 a nd Mct2 analysis for 
the SUSY(spsla) event when a pair of right handed 
squarks(mg Ji = 521 GeV) are produced, and each decays 
to a quark and an LSP(m x = 98GcV),(pp — > qaqR — > 
2x(q+Xi))- It is a parton level simulation with no initial 
state radiation(ISR) effect so that St ~ 0. 

Fig. [^a) is a scatter plot of the events in the My 2 - 
Mqt2 plane for x = m x- ^ s expected from the previous 
Mct a nd Mt analysis, most of the large Mt2 events 
are projected into a narrow endpoint region of the Mqt2 
variable, reflecting the large J ma x factor for the endpoint 
region of Mqt2 distribution. Fig. [2jb) is a compar- 
ison plot between the Mt2(x) anc > Mct2(x) distribu- 
tions for x — l-5m x [green], 2m x [bluc], 2.5m x [rcd], and 
3m x [yellow] . The uprisen distributions are Mqt2 (x) ~ X 
and the laid distributions are Mt2{x)~ X- F° r X = (1-5 ~ 
3)m x , the endpoint enhancement factors, J m ax(x) °f ©i 
are estimated to be 69, 33, 20, and 13. As in Fig. (TJ we 
can also observe a similar behavior such that if x is small 
then the Mqt2(x) distribution collapses into a region 
near zero with a diverging J maX y while it approaches to 
Mt2(x) when x na s a large value. 

As a result of a large J m ax factor, a small slope dif- 
ference in the Mt2(x) distribution can be amplified also 
in the Mct2(x) projection by a factor of J^ax; natu- 
rally transformed into a remarkably enhanced or uprisen 
endpoint structure of the mother particles. Therefore, as 
discussed in the previous section, the Mct2(x) projection 
also can help us to extract a meaningful endpoint defined 
as a BP in the Mt 2 (x) distribution. Even though there 
exist heavy systematic uncertainties in the signal and 
backgrounds, the dim BPs will be amplified and it can 
significantly reduce the systematic uncertainties we suf- 
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FIG. 2: (a) M T2 vs. M CT 2 of SUSY(spsla) q R (- 
pair production with a trial WIMP mass of 100 GeV, and 
(b) Mct2(x) ~ X (uprisen) and Mt2(x) ~ X distributions for 
X = 1.5m x [green], 2m x [blue], 2.5m x [red], and 3m x [yellow]. 



fcrcd from in fitting with more elaborated fitting schemes. 
Fig. [3] shows the Mqt2 (uprisen) and M^laid) distribu- 
tions of slepton pair production events, where a pair of 
right(lcft) handed sleptons(m; H ^ = 191 (256) GeV) arc 

produccd(sps5), and each slepton decays to / + Xo with 
m^i =119 GeV. This is also a parton level result and 
we used \ ~ m x},' I R the Mt2 distribution, there exists 
a BP which is made by the endpoint of Ir but buried in 
the II signature. One can see that such a BP is enhanced 
in the Mqti projection, where the J m ax (x) factor for In 



FIG. 3: The amplified endpoint structure of the RH and 
LH sleptons in Mct2 (uprisen) distribution and usual Mt2 (x) 
distribution(laid) 



is 5.3 for x = m x\- With this example, one might natu- 
rally think about the possibility to resolve some mass hi- 
erarchies with reduced systematic uncertainties using the 
Mct2{x) projection of some inclusive event data. The 
mass hierarchies will be represented by several po values, 
for instance, here there are two po values to be resolved 
as follows 



2mr 



(16) 



each of which can be extracted from the first and the 
2nd endpoints of the Mct2(x) or Mt2(x) distributions, 
respectively. 

As we will sec in the next section, Mct2{x) projec- 
tion can be a very powerful tool for resolving several 
| Po| values originated from various mass hierarchies be- 
tween new particles. Basically, wherever there exist sym- 
metric decay chains with various mass hierarchies, one 
can try to resolve them in terms of various |po| values. 
If the symmetric chain involves three body decays like 
(2 x (Y — > qqX)), the effective |po| value which con- 
tributes to the endpoint of Mct2(x) or Mt2{x) will be 
(my — mx)/2. (Here the Mct2(x) an d Mt2(x) is al so 
defined with only 2 visible particles with properly modi- 
fied missing Ifi T .) From now on, we define |po| by 



|po| = Maximal absolute momentum of daughter 
particles in the rest frame of mother particle 

(m Y ~ m|)/(2my) for 2 body decay 
(my — mx)/2 for 3 body decay. 



(17) 
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Therefore, the Mqt2(x) projection enables us to imple- 
ment some precision endpoint measurement using the in- 
clusive signal and background data with large uncertain- 
ties. The amplification of some BPs by Mct2(x) P ro ~ 
jection increases our observability for meaningful signal 
cndpoints in Mt2 distribution, which are buried in severe 
backgrounds. 



IV. APPLICATIONS 

In this section, we will describe more useful example of 
Mct2(x) projection in order to obtain some constraints 
between mg, m^, and other superparticle masses in a sit- 
uation where large uncertainties are involved with jets. 
Specifically, an example of measuring some superparti- 
cle mass constraints via Mct2 projection is shown, us- 
ing 6 hard jets + Ifi T signal from a SUSY model. The 
mSUGRA type SUSY benchmark spectrum is chosen in 
Table H 



N-jet + Missing E T events from qq production & 

important 2-jet sub-M T ,/M CT , systems 



ISR 



rriq 


m- g 


m.± ,-a 


m x? 


1026.3 ~ 1036.6 GeV 


649.4 GeV 


182.1/181.2 GeV 


98.6 GeV 



TABLE I: A benchmark SUSY spectrum. SOFTSUSY[27| is 
used to calculate the spectrum with mSUGRA model param- 
eters, my 2 = 250 GeV, M = 900 GeV,tan/3 = 10, fj, > 0, and 
A = 0. 

At this benchmark point, squark and slepton masses are 
heavier than gluino by a few hundreds GeV. Also the 
2nd ncutralino \2 an( i chargino xt are mostly wino so 
that their masses are nearly degenerate. Thus, if we take 
into account the production and decay rates at the LHC 
energy of 14 TeV, the important new particle mass hier- 
archies can be categorized by 4 mass scales listed in Table 
HI gg and gq productions are the dominant superparticle 
generation processes with a ~ 3 pb for both cases, qq 
production also has a sizable cross section about 0.5 pb. 
Lepton production is suppressed here because sleptons 
are also very heavy so that leptons cannot be generated 
via some cascade decay chain of gluino or squarks. So 
most of new physics signals are hard N-jets + lfi T signal 
from squark or gluino decays. Several important squark 
and gluino decay chains give the following branching ra- 
tios 

1. BR(q -> g + q) ~ 65% 

2. BR(g ^Xi+Q + q)- 44% 

3. BR(g -^x° 2 + q + q)~ 25% 

4. BR{g -^Xi+q + q)- 15% 

Here, all the listed gluino branching ratios are for 3 body 
decay processes via virtual squarks. The inclusive hard 
N(> 6) jets + Ifirp signal can come from heavy squark 
pair production where each squark decays to j s + g, and 



1037 q 



649 




IS2 



J : Decay mode 

Jet from a 
decay mode 

^ : Missing particle 
: 2-jet subsystem 
^ of interest ^ 

*W,)=314 
1 st endpoint 

w s .)= 275 

Sizable but small BR 
->P (/ JfJ ) = 234 
2° d endpoint 



FIG. 4: A SUSY event topology for hard N(^ 6)-jets + $ T 
signature and the important 2-jet Subsystem- Mt2 and Mcti 
systems (yellow boxes). j s : a jet from q decay to g, j g \ : a 
jet from g decays to Xi> 3g2 '■ a j et from g decay to X2/X1J 
and j x : a jet from X2/X1 decay to Xi- 



subsequently, g —> jgjg+Xi/xf /x.2 through 3 body decay 
chain, 

pp -> G(—6 T ) + q + q 

G(-5t) + 2 x ( 3s + ~g(^ j g j g xVxt/xl)), (18) 

xtlxl^w/z*(-*j x j x ) + x ° u 

where G(~8t) is the extra particles(e.g. ISR) which gives 
a recoil transverse momentum St to the qq system. Fig. 
|4] shows the event topology of the qq production we con- 
sider. j s and j g in (|18[) denote the jets from squark decay 
and gluino decay, respectively. j g can be separated to j g i 
and j g 2 as in Fig. 01 according to their decay processes. 
In this paper, xt an d X2 are assumed to decay to light 
jets(j' x ) + Xi via W or Z bosons if they are produced in 
the gluino decay. 

In this regard, phenomcnologically the most important 
question would be, "How can we measure the superpar- 
ticle masses using N-jets + Jfi T signal in a situation with 
large jet uncertainties in identifying true signal jets with 
complex combinatorics?" One might attempt to measure 
the gluino or squark masses using the so-called 'M^-kink 
method [H|'. According to this method, one can have a 
good chance to measure the squark(or gluino) and x\ 
masses simultaniously by extracting the position of the 
kink, if 6 jets(or 4 jets) from squark(or gluino) pair decays 
are efficiently identified. In particular, for gluino mass 
measurement, if squarks are very heavy so that they are 
decoupled at the LHC, then most of the hardest 4 jets 
would be from the decays from the gluino pair, the ef- 
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ficiency of choosing the correct 4 jets can be quite nice. 
However, when squarks are heavier than gluino by just a 
few hundred GeV, then the efficiency gets worse, partic- 



ulary for Pr(j s ) > Prijg) when 



rn g ) is comparable 



to, or larger than (nig — m^o). Although the squark pair 
production rate decreases as the (nig — m g ) increases, 
the two j s with hard Pt are more likely to be selected in 
the gluino Mti calculation, and badly pollute the end- 
point structure (T3 |. Furthermore, the ISR can be an 
important source of jet backgrounds with hard Pt also 
[HI, . Thus, there exists a large amount of systematic 
uncertainty in correctly estimating the jet backgrounds 
of gluino Mti ■ 

As we will describe soon, the Mqt2 projection can 
help resolving several mass differences hidden in various 
inclusive jet signatures in general. The mass differences 
are represented by several |po| values in the endpoints 
measurements of the Mti and Mqti distributions, as 
we discussed in the previous sections. Such |po| values 
are listed in Fig. FJ] whose values have been calculated 
by (TIT]). 

Here we concentrate on measuring the mass difference 
between squark and gluino since the usual gluino mass 
measurements have a lot of uncertainties as we explained 
in the previous paragraph. After the description of this 
specific example, we will present a generalization of the 
Mcti amplification for resolving several mass hierarchies 
which have been buried in various inclusive signatures, so 
regarded as meaningless in the Mti distribution. 

Our strategy for measuring the squark-gluino mass 
difference is naturally focusing on the two j s rather 
than trying to select four j g from gluino pair decays. 
Surely, there exists a large amount of uncertainty in cor- 
rectly choosing two j s and estimating the backgrounds 
also. However, we attempted to calculate Mct2(x) an d 
Mt2(x) f° r the Subsystem (qq — > j s jsgg), considering 
two gluinos as the effective missing particles for the sys- 
tem. Once we could construct the Subsystem- Mt2{x) 
(M. Burns et. al. in 14j) using the correct pair of two 
j s , then in the limit of vanishing St, its endpoint corre- 
sponds to M 2 in © with |p | given by 



\Po(j s )\ 



2rriq 



314.2 GeV. 



(19) 



Therefore, the measurement of the Subsystem- Mt2(x) 
endpoint can provide us a constraint between nig and 
m g . Once we assume that the nig and m^o masses could 
be already obtained through the Mr 2 -kmk methods us- 
ing 6 hardest Pt jets, then m g might be also measured by 
(|19l) . However, as commented previously, the endpoint 
measurement of the Subsystem- Mt2(x) distribution is 
not easy. Choosing the correct pair of two j s suffers from 
the jet uncertainties, and hence some effective event 
selection cuts or methods purifying the signals should 
be applied. Thus, we describe another approach to the 
endpoint measurement of the Subsystem-My 2 (x)- In 
short, it can be summarized as follows : 



We do not care whatever backgrounds there exist. 
Once there is a dim BP from a signal endpoint in the 
Mt2 distribution with large systematic uncertainties, 
then we can try to extract the position of the amplified 
BP in the Mqt2 projection. 

With this point of view, we calculate the Subsystem- 
Mct2(x) an d Subsystem- Mt 2 (x) variables for all pos- 
sible 15 pairs of the two jets among the 6 highest Pt 
jets in an event with N(= 6,7)-jets + Jfi T . If the event 
is from a real squark pair production, then there exists 
at least one true pair of two j s among the 15 pairs of 
possible two jets, and it might consistently contribute to 
the meaningful endpoint in the Subsystem- Mt2 (x) dis- 
tribution, although the slope change must be faint. The 
portion of true two j s pair will be much more reduced if 
we take into account the background events, e.g. gg, gq 
+ etc with hard ISR effects. Anyway, using the events 
with N(= 6, 7)-jets + Jp T , we tested the possibility to 
measure the corresponding amplified BP in the Mct2 
projection with less systematic errors. Extracting mean- 
ingful endpoints using more inclusive events with general 
N(> 6)-jcts + Ifirp will be discussed in the latter part 
of this section and a forthcoming paper [25[ . We found 
that even with this small fraction of N(= 6, 7)-jets + Ifi T 
events, a sizable number of events survive the usual new 
physics cuts, and contribute to some meaningful ampli- 
fied BPs of the inclusive Subsystem(IS)-McT2(x)- To 
reconstruct the BPs in the IS-Mct2{x) an d Mt2(x), the 
15 values of Subsystem- Mct2{x) an d Subsystem- M^ 2 (x) 
are constructed with given trial gluino masses x f° r an 
event as follows and we made histograms of all the calcu- 
lated values without any preferential jet selection among 
the 6 hardest jets, 



M, 



(n) 

CT2,T2 



m£t%{xY 



(x) 



min [masi{M^ T 1 l,M^f T }](20) 



X - 



(n)2 
(n) 



2(e 



(«), 



±P$-k«0 



The index n(= 1, ■ ■ ■ , 15) means a specific combination 
to select two j s among the 6 hardest Pt jets. p\ T ' are the 
transverse momenta of the two jets selected as j s in the 
lab. frame. J2 T is the total missing transverse momenta, 

and is modified missing transverse momenta for 

the subsystem with indicating the pt sum of the 
remaining 4 jets which are selected as four j g in the n-th 
combination. \ 1S meant to be a trial gluino mass with 
this setup of visible and missing momenta. 



We simulated the qq, gg, qg production events of 
10 fb" 1 at the LHC energy of 14 TcV using PYTHIA 
6.4 [26| with ISR/FSR turned on. Fully showered and 
hadronized events were passed to the PGS 4.0 detector 
simulator |2n . The energy resolution parameter in the 
hadronic calorimeter was given by AE/E = 0.6/VE, and 
jets were defined using a cone algorithm with AR = 0.5. 
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x/|Po(>)l 


Jmax (x) 




M^ 2 ax (xT xp 


a ,rmax i \true 
M CT 2 (X) 


71 /rmax i \ exp 
M C T2 (X) 


gsys /g s y s 




1.24 
2.57 


12.2 
3.9 


814.8 GeV 
1179.8 GeV 


810.2 ±19.5 ±31.6 
1150.2 ± 18.3 ± 30.4 


518.6 GeV 
998.5 GeV 


520.5 ± 1.6 ± 0.3 
1003.1 ±4.5 ± 0.8 


0.01 
0.03 


0.12 
0.11 



TABLE II: Expected/Fitted endpoints of the Mct2 and Mt2 for x — l-24|po(j s )|, 2.57|po(j s )|- For each measured endpoint 
value, the first error is statistical and the second one is systematic in the fitting. One can find that statistical errors are reduced 
by a factor of 1/ J m ax(x) as expected, while the the systematic errors are reduced by a factor of 1/ Jmax{x)i or l ess - 



To be more reliable for the new physics event experi- 
ments, we imposed several cuts for pure N-jets ± Tfi T 
events [l], 0, H^l although we did not include the SM back- 
ground event samples for simplicity. 

The event selection cuts were as follows: 



1. No leptons, no b jets in the event, 

Number of jets = 
P T th > 50 GeV, 



2. Number of jets = 6, 7 with P T st > 2nd > 100 GeV, 



3. # T > 100 GeV, 

4. a n > 0.45 with n(= 1, 
of selected two jets, 

5. A T (= |# T ±V, 



15) for the pairs 



, 6 PtI 



< 30 GeV, 



where the index- j denotes the 6 hardest jets. The 
a n = P^/ 171 ]] UM where P^ 2 and m ™j are the Pt of 
the second hardest jet and dijet invariant mass, respec- 
tively, for the n(= 1, • • • , 15)-th pair of selected two jets. 
If gg decay directly to XiXi> then At becomes St which 
is the total transverse momentum of the squark pair sys- 
tem defined in (fT5|) . As explained in the previous sec- 
tion, the fifth cut is imposed to suppress the shift of 
the Mct2{x) endpoint (|15[) under arbitrary transverse 
boost of two mother particle systemic?). 4 We calculated 

the Mpy 2 (x) and M T n 2 (x) f° r &H the events passing the 
cuts. The single jet invariant masses, in calculat- 



and A/py 2 °f (EHJ) wa s ignored since we 



ing the M$ 

found that it is quite helpful for the endpoints of both 
variables to be located at the expected position, reducing 
the jet energy resolution effect. Ignoring the jet masses 
also satisfies all the inequalities of Mt2 and Mct2 and 
the constructed variables can saturate to the boundary 
value because of enough statistics of light QCD jets in 
the events. In our SUSY event sample, the jet multiplic- 



ity ratio is 



r(jV Jet =6,7) 



0.23. 



With the relatively small 
= 6,7) surviving the 



a(A^ t >6) 

portion of the event sample (Nj et 
cuts, we still could observe the BP in the IS-Mt2(x), and 
the amplified BP in the IS-Mct2(x) also. 

Fig. (5Ja,b)(two histograms in the top) shows the 
Mt2 and the corresponding Mqt2 distributions for x = 
1.24|p (j s )| = 389.7 GeV. Fig. |3c,d)(two histograms in 



4 The Mct2 endpoint with a sizable St will be studied in a forth- 
coming paper 12511 . 



the bottom) are also Mt2 and Mct2 for x = 2.57 |po| = 
806.5 GeV. The endpoint enhancement factor, J mQX 
is 12.2 for (a)-Kb), and 3.9 for (c)^(d). We con- 
structed many histograms with respect to various \ val- 
ues. Among them, we selected the two cases as shown 
in Fig. [5] Actually, compared to |po|, a much larger or 
smaller value of x is n °t a good choice as explained in the 
previous section. Based on our trial experience, choosing 
X of order of |po| shows a proper amplification feature 
with clean breakpoint structures, and this selection can 
be realized by observing an average value of \Pt\ of the 
hardest jets. A strong A^-supprcssion cut was quite ef- 
fective for observing a BP with a small slope discontinuity 
at the expected endpoints in the Mt2 distributions, and 
we could observe the amplified BPs also at the expected 
position in the corresponding Mct2 distributions. Actu- 
ally, this cut is effective for selecting the qq production 
event. The expected endpoints are pinpointed by a red 
dashed line in each figure. The bin interval was selected 
as the best one among several bin candidates, by which 
the histogram near the expected BP regions shows its 
characteristic feature maximally while keeping its statis- 
tical relevance. 

In Fig. EJa), the Mt2(x — 389.7) distribution shows 
a small slope discontinuity near the expected endpoint 
814.8 GeV. One can find that the faint BP structure is 
amplified near the expected endpoint 518.6 GeV in the 
Mct2(x = 389.7) (Fig. [5jb)) projection, producing a 
sharp cliff wall. Similarly, there exists a small BP near 
the expected endpoint 1179.8 GeV of M T2 {x = 806.5) 
distribution in Fig[5tc), and the BP is transformed to 
another BP with a sharp cliff near the expected endpoint 
998.5 GeV of Af CT2 (x = 806.5) in Fig. [S^d). 

Table |TT] shows the expected and fitted endpoints of 
the Mt2 and Mqt2 distributions for x = l-24|po | and 
2.57|po|. The errors after the fitting process are also 
listed in the columns of M T n £ x (x) exp and Mg^f (x) exp - 
Here, the first error represents statistical one, and the 
second represents the systematic one in the fitting pro- 
cess. The fit model functions for the endpoints were a 
Gaussian smeared linear function for a signal, and a lin- 
ear function for the backgrounds in the Mt2 distributions 
(Fig. ETa,c)), 



/(to) = 9(to- M^ x ) 



2na 



T2 Q,- m ')2 , 

e ( — > 



x ai(m' - MT2 X )dm / + a 2 m ± a 3 . 



(21) 



Also, the endpoint region of the Mqt2 distribution (Fig. 
[5£b,d)) was fitted by the Gaussian smeared step functions 
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for signals, and step function for backgrounds, 

„, , ai0(m - M rT9 y ') r ni cT2 (m _„ t ; 
/(m) = — i — _£T2 — ; y e ( — ^2—) dm ' 



27TCT 



a20(m — M CT2 



27TCT 



i CT2 (m-m") 2 * 

e ( 5^ — )d m " 



(22) 

Fitting was implemented by the binned x 2 method us- 
ing MINUIT with MINOS error processing [HI which 
takes into account both parameter correlations and non- 
linearities. The two step functions for signal endpoints 
was designed in (|22j) to fit two cliffs appearing in the 
Mct2 distributions. We will explain the second end- 
point below the 1st endpoint which were pointed by red 
dashed line in the Mqt2 distribution (Fig. E][b,d)) in the 
latter part of this Section. For simplicity, we used the 
same smearing widths a for both the two endpoint cliffs 
in the Mqti distribution. Actually, the uncertainty from 
the endpoint broadening effect by some widths is not our 
issue here. As described in Sec. HH it is simply because 
the endpoint uncertainty reduction factor is expected to 
be just the order of 1/ J max in that case, and the reduced 
width is compensated by the explicit error propagation 
factor J max for obtaining |po| value. It is natural that 
there is no advantage in using Mcti projection for re- 
ducing the statistical errors. Anyway, the fitted endpoint 
values were found not to be very sensitive to the choice 
of a in the range 0-10 GcV for M T2 , and 0-1 GeV for 
McT2, and we fixed the smearing widths a to be 5 GeV 
for the Mt2 endpoint, and 0.5 GeV for the Mqt2 end- 
points. The width for the Mqt2 endpoint is a suppressed 
value from the width of the Mt2 endpoint by 0(1/ Jmax)- 

The green lines in Fig. [5] represent the fitted model 
function. In order to estimate the systematic uncertainty, 
we attempted to fit various endpoint regions with the 
simple model functions while keeping the % 2 /ndf < 2 af- 
ter fitting( ndf = the number of fitted bins - 1). The 
small value of x 2 /ndf might guarantees the plausibility 
of the fit range described by the model function. Ac- 
tually, for anyone trying to find some optimal range of 
the fitting in the Mt2 distribution, he easily notices that 
there exist a large number of choices for a lower boundary 
and an upper boundary, which is consistent with small 
X 2 /ndf after fitting. In this situation, because the BP 
is quite a faint one can observe that the fitted endpoint 
parameter changes very much with a large uncertainty 
as the fit range of validity is varied. On the other hand, 
the fitting of the Mqt2 endpoint regions can be different. 
Due to the emergence of the sharp cliffs, which are the 
amplified BPs, we could clearly elaborate our effective 
fit model functions and the fit ranges of validity. Ac- 
tually, there is another meaningful BP near 690 GeV in 
the Mt2(x = l-24)po|) distribution, which is now clearly 
seen as the cliff near 505 GcV in the Mct2 projection(For 
X = 2.57|p0|, a similar BP enhancement is observed). It 
looks even more clearer than the first BP we expected 



for the squark pair decays. In the past, one would have 
neglected such a BP buried in the region of the low Mt2 
distribution, as well as suffer from large systematic un- 
certainties in extracting the position of the first BP. Now, 
in addition to the first cndpoint(red dashed), we can try 
to extract the amplified second lower endpoint and inter- 
pret its origin with less uncertainty. 

One can easily see that the second endpoint is originat- 
ing from the mass differences between the g and X2, Xii 
represented by another |po| value, 



|Po(j 3 2)| 



234 GcV. 



(23) 



This is the maximal momentum of jets from a gluino 
three body decays to j g j g + X2/X1 m the rest frame of 
the gluino. Since the BR(g —> j g 2jg2+X,2/xf) is 4-6 times 
larger than the BR(g — > j g ij g i + Xi)j gluino decays to 
X2/X1 can contribute significantly to the inclusive N(> 
6)-jets + Jfirp signature than the gluino decays to x? since 
most of the X2/X1 hadronically decays to 2-jets(j x ) + Xi 
via Z* /W bosons. The point we want to emphasize here 
is that, once there is some possibility of a pair of mother 
particles to be symmetric in the lab. frame(^ Pt = 0), 
then as analyzed in Sec. [TTJ it might contribute to some 
other endpoint in the M^ 2 (x) distribution with some |p | 
value defined in (jTTJ) . For example, in the squark pair 
production events we consider, the symmetric condition 
of a pair of gluinos can be realized if 



(24) 



1 effective ^? T , if the two p£p in (f20|) correspond to 



However, regardless of the gluino pair symmetry condi- 
tion, M$ and M^ 2 in d^HD 

can include Subsystem- Mt2 
and McT2 for the subsystem which consist of two visible 

.(n 
Jg2 

the transverse momenta of two . Here means the 
2 gluino jets j g 2, each of which comes from each mother 

(n) 

gluino, and is selected for p iT . The corresponding sub- 
system is described in the third yellow box of Fig. [4] If 
denotes the other two jg2 which are not selected as 
the subsystem visible particles, then the effective miss- 
ing transverse momenta, JS^ in the (|2"0|) , and At in the 
event selection cuts become as follows: 

= %, 



<rp ~r Q J- 

{TPtU 



E^t(j x )-M 



+{E / m^)+E / m42 ) )} 

= - E ^(42 ) - E P ^x) + St, (25) 
A T = I#t + E j=1 .. 6 PtI 

= \5t-J2 P tUx)1 (26) 

where the assumption that j x is not included in the 6 
hardest Pt jet candidates is adopted. Then, our strong 
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FIG. 5: Endpoint fittings of the Mti (left columns (a) and (c)) and Mct2 (right columns (b) and (d)) distributions for 
X = 1.24|po(j' 3 )| ((a) and (b), top row), and for \ = 2.57|po(Js)| ((c) and (d), bottom row). The green lines are the fitted 
model functions and the red dashed lines indicate the expected endpoints : 

M£ir(1.24|p |) = 814.8 GeV, Mg?£ (1.24|p |) = 518.6 GeV, M^ 2 aa: (2.57|po|) = 1179.8 GeV, M$?f (2.57|p |) = 998.5 GeV. 



Ay suppression cut implies that the effective missing 
transverse momenta are 

^«)~-£p T (4«)). (27) 

This indicates that Subsystem- Mti and Mct2 m ([20)) 
can have a good chance of contributing to the well-defined 
endpoints described by | Po ( Jg2 ) I value if the gluino sym- 
metric condition (f2~4"|) is satisfied. In that case, the con- 



tribution consistently forms the slight slope change in 
the Mti distribution, and it will emerge as a sharp cliff 
in the Mqt2 distribution. In our example, the strong 
Ay suppression cuts of (|2l)1) effectively acts as a St sup- 
pression cut, and then it will create a bias to choose 
small ^PtUs) events because ^PtUs) distribution is 
expected to be an even convex function when St is zero. 
Then, the symmetric condition can be effectively satisfied 
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A r (JV) = 5 T : Total transverse- 
momentum of N-th shell 



A < iV - 1) 



: Tr. momentum of (N - l)-th shell 



A r (N) ~ 2j PtUi (N ' <w= 1 '" 2,vC2) ) 
;=i,2 



T,T, : m-th mother particle pair, 
symmetrizable if A T ' m) = 0, m = l..N 

• : 2-jet(/= 1,2) candidates 

for ra-th shell with the index i m = 1 ■ - 2m C 2 , 
subsequently defined from m — N. 



set of jets g {/"' ! "- , |w<n<7V} 



In m-th shell, 

M^ 2 Cr)-min [ max {M ( ^/"\ M 

A=6 



(2Xm,i M )\ -i 
CTIT i J 



A = k 1T + k 2T - X T <m,i '" ' with X T <m,im) = K +^P T G (m ''" ,y ) 
A T (m> = zf/"' +y - £ P T Ui (m+l,im ^) 

'=1,2 



FIG. 6: Extension of the Subsystem- Mct2,Mt2 in (|20[) for m-th(m — 1..N) symmetrizable mother pairs with symmetrization 
cut parameter A™ using inclusive hard n(> 2iV)-jet + lp T signatures. 



for a lot of events, and there exist large contributions to 
form some slope change near the expected exact second 
cndpoint. 

However, one can easily notice that the At suppres- 
sion cut cannot guarantee the exact symmetric condition 
of the gluino pair. In order to access to each symmet- 
ric mother particle pairs possibly hidden in the N-decay 
shells, one can extend the definitions of Subsystem- Mti 
and Mct2 of (|20|) as shown in Fig. [6] Here, with hard 
n(2N < n < 4N)-jet + Ifi T inclusive signatures, one can 
subsequently construct all the possible m-th Subsystem- 

M^™2)-^T2 variables from the highest shell, which has 
a well-defined generalized symmetric condition parame- 
ter. A™. In general, Subsystem- M^) 2 and Mj^ have 
C al ( = 2mC2) multiplicity which is the number of jet 
combinations possible in the m-th shell. As a result one 
can explicitly control those cut parameters to observe the 
amplified BPs from several symmetrized mother particle 
pairs. Furthermore, one can find that to realize the hid- 
den endpoints the one-dimensional decomposed version 
[H of Subsystem- M$ 2 and M^\ x and M^j, 

can be useful, if it is constructed using only partial com- 



ponents of PtU^) perpendicular to A^™ direction. In 
this regard, a more general Mqt2 spectroscopy will be 
discussed in a forthcoming paper [251 1 . As a result sev- 
eral mass hierarchies in various new physics models can 
be probed by general Subsystem- Mct2{x) amplification 
using the inclusive N(> 2, 4, 6, • • • )-jets/leptons /photons 
+ Tp T signatures, with less systematic uncertainties. 

Finally, we summarize the fit results of the first end- 
points(red dashed) in the column, M^ l 2 ax (x) exp and 
Mgffi (x) exp , of Table HI Using dozens of fit results 
for various fit ranges, we evaluated the statistical errors 
as 19.5 GeV for M T2 (x = 1.24|p |), and 1.6 GeV for 
Mqt2(x — 1-24 1 po|). The listed systematic uncertainties 
were obtained by la from the fitted endpoint values. It 
is noticeable that the systematic errors are suppressed in 
the Mgpf (x) exp by the factor of 1/J^ ax for X = l-24|p ] 
or less than that for x = 2.57|po|, while all of the sta- 
tistical errors are reduced by the order of l/J m ax- As a 
result, the precision to measure the |po| can be enhanced 
at least by a factor of J max as listed in the last column 
of Table |TTJ These results are well matched with our ex- 
pectation of the endpoint structure enhancement by the 
Mct2 projections. 
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Including more various systematic error effects (the SM 
backgrounds, the signal and BGs in higher orders, PDF) 
is beyond the scope of this paper. But, the main point of 
this paper is clear: The features of the BP enhancement 
by Mqti projection has a more stronger power in re- 
solving the buried mass hierarchies as our situation gets 
worse with large uncertainties. 

V. CONCLUSION 

We have introduced the constransverse mass, Mqt2, 
and have shown that it can significantly increase our 
ability to measure the endpoints precisely by amplifying 
the small slope discontinuity of the endpoint region, 
which are usually buried in the complex backgrounds 
with large uncertainties. As a result, the proposed Mqti 



variable for various subsystems with two visible particles 
can resolve the mass hierarchies hidden in various 
inclusive signatures for the event type, N(> 2, 4, 6, ■ ■ • )- 
jets/leptons/photons + Ifi T . We hope that this variable 
can be a useful tool in the LHC era. 
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